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ABSTRACT 



The quantitative spectral analysis of low resolution (~ 5 A) Keck LRIS spectra of 
blue supergiants in the disk of the giant spiral galaxy M81 is used to determine stellar 
effective temperatures, gravities, metallicities, luminosites, interstellar reddening and 
a new distance using the Flux-weighted Gravity-Luminosity Relationship (FGLR). 
Substantial reddening and extinction is found with E(B-V) ranging between 0.13 to 
0.38 mag and an average value of 0.26 mag. The distance modulus obtained after in- 
dividual reddening corrections is 27.7 ±0.1 mag. The result is discussed with regard 
to recently measured TRGB and Cepheid distances. The metallicities (based on ele- 
ments such as iron, titanium, magnesium) are supersolar 0.2 dex) in the inner disk 
(R < 5 kpc) and slightly subsolar (ps -0.05 dex) in the outer disk (R > 10 kpc) with 
a shallow metallicity gradient of 0.034 dex kpc -1 . The comparison with published 
oxygen abundances of planetary nebulae and metallicities determined through fits of 
HST color-magnitude diagrams indicates a late metal enrichment and a flattening of 
the abundance gradient over the last 5 Gyrs. This might be the result of gas infall 
from metal rich satellite galaxies. Combining these M81 metallicities with published 
blue supergiant abundance studies in the Local Group and the Sculptor Group a galaxy 
mass metallicity-relationship based solely on stellar spectroscopic studies is presented 
and compared with recent studies of SDSS star forming galaxies. 

Subject headings: galaxies: distances and redshifts — galaxies: individual (M81) — stars: 
abundances — stars: early-type — supergiants 
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1. Introduction 



The determination of the chemical composition and distances of galaxies is crucial for 
constraining the theory of galaxy formation and evolution in a dark energy and cold dark 
matter dominated universe. Ultimately, these measurements lead to ever stronger constraints 
on the cosmological parameters and the history of cosmic chemical enrichment, from the 
primordial metal-free universe to the present-day chemically diversified structure. For instance, 
the relationship between central metallicity and galac tic mass appears to be a Rosetta stone to 



understand chemica 



Maiolino et al. 



evolution and galaxy formation (|Lequeux et al. 



1979 



Tremonti et al. 



2004; 



2008). In a similar way, the observed metallicity gradients in s piral galaxies. 



apparently large for spirals of lowe r mass and shallow for high mass galaxies (IGarnett et al 



1997 



Skillman 



1998 



Garnet 



2004), provide crucial insight into galaxy formation and evolution. 



Both the observed mass-metallicity relationship and the abundance gradients are used to test 
the theoretical predictions of hierarchical clustering, galaxy formation, merging, infall, galactic 
winds and variability of star formation activity and IMF obtained in the framework of a 



ACDM dominated universe (Prantzos & Boissier 



2008 



Yin et al 



2009; 



Sanchez-Blazquez et al 



2000; 



* 



Naab & Ostrikei 



2009; De Lucia et al.1 



2006; 



2004; 



Colavitti et al 



de Rossi et al 



2007; 



Finlator & Dave 



2008 



Brooks et al. 



2007 



Koppen et al 



2007 



Wiersma et al. 



2009; 



Dave et al. 



201 lailbh . Note that this is only a small selection of papers relevant to the subject, others are found 



in the references therein. 

However, as intriguing the observations of the mass-metallicity relationship and the 
metallicity gradients of galaxies are, the published results are highly uncertain. They rely on 
observations of H II region emission lines, mostly restricted to oxygen, and the analysis method 
applied is the so-called "strong-line method", which uses the fluxes of the strongest forbidden 
lines of (most commonly) [O II] and [O III] relative to Hp . Unfortunately, abundances obtained 
with the strong-line method depend heavily on the calibration used. As a striking example, 
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Kewley & Ellisonl (|2008f) have demonstrated that the quantitative shape of the mass-metallicity 
relationship of galaxies can change fr om very steep t o alm ost flat depending on the c alibration 



used. In the same way, as shown by 



Kudritzki et al. 



(2008() and 



Bresolin et al 



(2009) in their 



study of the Sculptor spiral galaxy NGC 300, metallicity gradients of spiral galaxies can change 
from steep to flat and absolute values of metallicity can shift by as much as 0.6 dex, again as the 
result of different calibrations of the strong line method. In consequence, galaxy metallicities 
are uncertain by 0.6 to 0.8 dex because of the systematic uncertainties inherent in the strong line 
methods used. This major problem requires a fresh approach and is begging for the development 
of a new and independent method less affected by systematic uncertainties. 

An obvious alternative method to constrain metallicity is the detailed quantitative 
spectroscopic analysis of individual blue supergiant stars (BSGs) in galaxies. BSGs of spectral 
type A and B are massive stars in the mass range between 12 to 40 M in the short-lived 
evolutionary phase (10 3 to 10 5 years) when they leave the hydrogen main sequence and cross 
the HR-diagram at constant luminosity and almost constant mass to become red supergiants. 
Because of Wien's law massive stars increase their brightness in visual light dramatically when 
evolving towards l ower temperat ures and reach absolute visual magnitudes up to My ~ -9.5 mag 



in the BSG phase (IBresolin 



2003), rivaling with the integrated light of globular clusters and dwarf 
galaxies. Because of their extreme brightness they are ideal tools to accurately determine the 
chemical composition of young stellar populations in galaxies. 

BSG spectra are rich in metal absorption lines from several elements (C, N, O, Mg, Al, S, Si, 
Ti, Fe, among others). As young objects with ages of 10 Myrs they provide important probes of 
the current composition of the interstellar medium. Based on detailed high resolution, very high 
(S/N) studies of blue supergiants , which yield abundanc e s as accurate as 0.05 d ex 



signal-to-noise 


(S/N) 


(Przvbilla et al. 


2006 



Schiller & Przvbilla 



2008 



Przvbilla et al.ll2008ah . iKudritzki et al. 



;2008) 



developed an efficient new spectral diagnostic technique for low resolution spectra (FWHM ~ 5 
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A) with good S/N ratio (50 or better), which allows for an accurate determination of effective 
temperature, gravity, metallicity, interstellar reddening and extinction. Metallicities accurate to 0.1 
to 0.2 dex for each individual target can be obtained at this lower resolution and S/N. The method 



has been applied to i rregular and spi ral galaxies in th e Local Gro up (WLM - 



Urbane 



Uetal 



a et al 



2008 



NGC 3109 



Evans et al 



2009) and beyond (NGC 300 - 



2007 



Kudritzki et al. 



C 1613 



Bresolin et al. 



Bresolin et al. 



2007 



2006; 



2008) 



; M33 - 



In this paper we present the spectral analyis of low resolution Keck LRIS spectra of 26 BSGs 
in the dis k of the giant spiral gal axy M81. M81 is one of the most massive spirals in the Local 



Volume (IMcCommas et al. 



2009). It has low foreground extinction with a galactic luminosity 



of 2.5 L* (corresponding to M^r = - 24 mag a nd = -23 mag) a nd is characteristic of disk 



galaxies seen at redshift surveys out to z ~ 1 (IWilliams et al. 



2009J). The star formation history 



and chemical evolution of this galaxy have been subject to extensive recent photometric studi es 



(Dal canton et al 



2009; 



Davidge 



2009; 



Williams et al. 


2009; 


Barker et al. 


2009; 


Durrell et al. 



2010). H II regions and Planeta ry Nebulae have been st udie d by 



1987) and 



Stanghel lini et al 



Stauffer & Bothur 



(2010) 



J( 1984L 



With 



extending the classical work by iGarnett & Shields! 1 
our work we provide for the first time direct quantitative spectroscopic information about stellar 
metallicity of the young disk population. 

An important additional aspect of the quantitative spectroscopy of BSGs is their use as 
accurate distance indicators through the Flux- weighted Gravity-L uminosity Relationship (FGLR). 



This new distance de termination method has been introduced by 



Kudritzki et al. 



Kudritzki et al 



(l2003h and 



(2008). It uses stellar gravity and effective temperature as a measure of absolute 
bolometric magnitude and provides a distance estimate which is free of the uncertainties caused 
by interstellar reddening, since the determination of reddening is a by-product of the quantitative 
spectral analys i s. Fir st distance d eterminations using this method have been carried out by 



Urbaneja et al 



(2008, WLM) and 



U et al. 



(120091 M33). 
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There has been a lo ng history of attempts to measure the distance to M81 from Hubble 



(1929) to the present (see 



McCommas et al 



.2009 for references and a plot of distance modulus 



as a function of time). The work published over the last decade gives a range between 27.60 to 
28.03 mag in distance modulus indicating an uncertainty of 20%. Our BSG spectroscopy and the 
effective temperatures and gravities determined will give us a FGLR distance which we can then 
compare with most recent HST work on Cepheids and the tip of the red giant branch, TRGB. M81 
has been used as one of the calib ration galaxies for th e Tully- Fisher and the su r face fluctuat ion 



methods in the HST Key project (|Freedman et al 



200 II) and by 



Mould & Sakail (120081 



2009) 



In section 2 of this paper we describe the observations and data reduction. Section 3 
discusses the quantitative spectroscopic analysis and the determination of extinction, effective 
temperature, gravity and metallicity. Section 4 and 5 discuss interstellar reddening and compare 
the spectroscopically determined stellar parameters with evolutionary tracks in order to constrain 
the evolutionary status of the objects observed. Section 6 compares metallicity and metallicity 
gradient of the BSGs with published metallicity constraints for the older disk population of 
M81 and discusses chemical evolution over the last Gyrs. In section 7 we provide a galaxy 
mass -metallicity relationship based on BSG spectroscopic studies and compare with published 
work using H II region emission lines. In section 8 we determine a new distance to M81 using the 
FGLR-method and discuss recent Cepheid and TRGB work. Section 9 sumarizes the results and 
discusses aspects of future work. 



2. Observations and Data Reduction 



The observations were carried out w ith the Keck 1 te lescope on Mauna Kea and the Low 



Resolution Imaging Spectrograph (LRIS, 



Okeet al. 



19951) using the atmospheric dispersion 



corrector, a slit width of 1.2 arcseconds, the D560 dichroic and the 600/4000 grism (0.63 A pix ! ) 
and the 900/5500 grating (0.53 A pix ! ) in the blue and red channel, respectively. In this paper, 
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we will discuss and analyze the blue channel (LRIS-B) spectra only, which have a resolution of 
5 A FHWM. Because of the UV sensitivity of the LRIS-B configuration the spectra extend to 
shortward of the Balmer discontinuity at 3640 A, which is crucial for the determination of T e ff 
from the Balmer jump (see section 3). Three MOS fields were prepared with 20 to 25 targets each. 
The BS G candidate targets w ere selected from HST ACS B,V images obtained within the ANGST 

I || I 

project (IDalcanton et al]|2009|) . which covers the whole galaxy. Published B,V photometry of the 
M81 ANGST fields was used to preselect targets with point source PSF characteristic and with 
-0.2 mag < B-V < 0.4 mag and V < 21.5 mag. Each target was carefully inspected with regard 
to multiplicity. Fig. Q] shows the selection from the color-magnitude diagram (CMD) and the 
location of our targets within the gal axy . None of our targets is relate d to one of the stellar clusters 



investigated by 



Chandar et al 



(1200 II) or 



Santiago-Cortes et al. 



QOlOj). 



The observations were scheduled for three dark nights in 2010 (February 14 to 16). The first 
night had perfect conditions with 0.75 arcsec seeing yielding reasonably exposed spectra with a 
total exposure time of 6.75 hours (observed in exposure segments of 45 minutes each) of the first 
field (field Z). The observing conditions degraded significantly during the second and third nights 
with poor seeing (1.3 arcsec) and occasional clouds. As a result, almost one half of these two 
nights was lost and only one additional field (field C) could be observed with a total of 1 1 .3 hours 
exposure time under mediocre conditions. 

Data reduction was performed using a custom pipeline written in IDL designed to efficiently 
extract faint objects observed over a full night. LRIS science and calibration frames were 
flat fielded and bias subtracted. For each reduced frame, object spect ra were trace d along 



the dispersion axis and extracted using the optimal extraction method (IHorne 



1986) meant to 



maximize the S/N of faint spectra. For this technique we utilized a Moffat function which 
was determined to best fit the 2-D spectral profile at each pixel (wavelength) perpendicular to 
the dispersion. The Moffat fit was modified to include a measure of the background level for 
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subtraction. The spectra were then wavelength calibrated using techniques in the idlspec2d IDL 
package developed for SDSS. 

Each science object spectrum was flux calibrated by performing corrections for wavelength 



dependent extinction at varying airmass over Mauna Kea (IBeland et al. 



1988) and then 



multiplying by a sensitivity function to convert extracted data numbers into units of ergs/s/cm 2 /A. 
The sensitivity function was calculated by scaling airmass-corrected observed flux standard s tars 



(GD 5 0, Feige 34, HZ 44, and BD+33d2642) to the published spectral energy distributions of 



Oke 



1990). A final spectrum for each target was produced by taking the median of all wavelength and 
flux calibrated spectral frames. Those spectra were normalized by manually selecting continuum 
regions and dividing by a high order polynomial fit to the continuum flux levels. The S/N values 
of our spectra vary between 40 to 80. 

Table 1 provides the information about the objects used for this spectroscopic study. While 
we selected 25 targets in each field, we could not use all of them. A few turned out to be blue 
foreground objects in the Milky Way halo, some had composite spectra indicating the presence 
of several objects in the slit and for some the S/N was not sufficient. For the remaining objects 
we list coordinates, galactocentric distance, spectral type, V magnitude, B-V color and the 
measured B aimer jum p Dg in Table 1. The way, how D# is defined and measured, is described in 



Kudritzki et al 



(2008). 



3. Spectroscopic Analysis 



The analysis method has been described in detail in 



Kudritzki et al. 



(|2008|) . A comprehensive 



grid of line-blanketed model atmospheres and very detailed NLTE line formation calculations is 



used to calculate spectral energy distributions (SEDs) 



synthetic spectra. Relative to the work presented in 



, including the 



3 aimer jump, and normalized 



Kudritzki et al. 



(|2008|) the grid has been 
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extended to cover temperatures from 16000K down to 7900K at gravities between log g =3.0 to 
0.8 (cgs). The lower limit of log g is a function of T e ff parallel to the Eddington-limit. Models 
are calculated for 14 metallicities [Z] = log(Z/Z ) : -1.30, -1.15, -1.00, -0.85, -0.70, -0.60, 
-0.50, -0.40, -0.30, -0.15, 0.00, 0.15, 0.30, 0.50 dex. Z/Z is the metallicity relative to the sun 

in the sense that the abundance for each element is scaled by the same fa ctor relative to its solar 

I — 1 I 1 

abundance. Solar abundances were taken from iGrevesse & Sauvall (11998b . except for oxygen 



where we adopt the value 
grid we refer the reader to 



rom 



Allende Prieto et al. 



Kudritzki et al 



(2001). For all further details of the model 



(|2008|) . The phys ics of the mode 



Przvbilla et al 



atmo spheres and the 



(2006) and references 



NLTE line formation calculations are described in detail by 
therein. 

The spectral analysis proceeds in several steps. First, fit curves in the (log g, T e ff)-plane 
are constructed, along which the models reproduce the observed Balmer jump and the Balmer 
lines. The Balmer jump is mostly a function of temperature, but also depends weakly on gravity, 
whereas the Balmer lines depend mostly on gravity and weakly on temperature. Fig. [2] and [3] 
demonstrate the dependence of the Balmer jump on temperature and of the Balmer lines on 
gravity. The intersection area of these fit curves determines the stellar effective temperatures 
and gravities and the corresponding uncertainties (see Fig. SI). The fact that the fit curves for the 
Balmer jump and the Balmer lines are not orthogonal leads to relatively large error boxes, in 
particular with regard to gravity log g. On the other hand, the flux weighted gravity 



log gf = log g-41og(r eff x 10 



-4> 



(1) 



is determined much more accurately, since the Balmer lines depend solely on log gp for 
temperatures higher t han 9000K (for an explanation of the physics behind this behaviour, see 



Kudritzki et al. 



2008). This is important for the use of flux weighted gravity as an indicator of 
absolute magnitude and distance (see section 8). Fig [51 [6j U\ M show fits of Dg and one Balmer 
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line for the remaining objects in Table 1 (with spectral types later than or equal to B3) to give an 
impression of the quality of the data. We note that we usually try to use all Balmer lines from H4 
to H10 to constrain gravity. However, varying from star to star we may encounter difficulties with 
individual Balmer lines. H7, for instance, is many times corrupted by interstellar Call absorption. 
H4, H5 and even H6 are sometimes affected by HII emission. Another problem are strong stellar 
winds, which can fill H4 and H5 with broad emission. Spectral flaws by improper corrections of 
comic ray hits may also affect line profiles. However, in general, we have more than one Balmer 
line per star to constrain gravity, usually three to four. For Fig. 2, 6, 8 and 9 we have selected the 
best fitting cases. 

Three objects of our sample are of earlier spectral type (BO. 5 to B1.5). For those, the Balmer 
jump is not a good temperature indicator. We use the ionizatio n equilibrium of Si II, S i III, and 



Urbaneja et al.l (|2005ah . which 



Si IV lines instead and apply the analysis method developed by 
relies on the use of line-blanketed NLTE model atmospheres including the effects of stellar winds. 
Fig |9] shows the spectral fits for the key lines of these objects. 

For three objects of later spectral type (Z4, Z9, C21) the wavelength range of the observed 
spectra does not cover the region of the Balmer jump. Thus, the only way to est imate their 



temp erature is the rel ationship betwee n effective temperature and spectral type (see lKudritzki et al. 



20031) . As shown by 



Kudritzki et al 



(|2008r) this method works only, as long as the metallicity 
is about solar. From the galactocentric distance of these objects and our study of metallicity 
and metallicity gradient (see section 6) for the other objects in our sample this seems to be a 
reasonable assumption and, thus, temperature, gravity and luminosity of these objects are very 
likely well determined. Nevertheless, we will not make use of these objects for the determination 
of the distance to M81 from the flux weighted gravity. 

We note that with the fit of either the Balmer jump Dg, or the silicon equilibrium, or the 
spectral type in the (log g, T e /y)-plane we can always calculate a reddening correction E(B-V) 
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along the fit curve at each effective temperature by comparing the observed value of B-V with 
the one calculated b y the atmospheric m odel. For fitting the Balmer jump, we then correct for 



reddening using the 



Cardelli et al 



(|1989|) reddening law with Ry =3.1. Once we have found the 
intersection with the fit curve for the Balmer lines, we can then also determine the final reddening 
value E(B-V) for the final values of T e ff an d l°g g- It is a big advantage of this spectroscopic 
determination of stellar parameters that it yields interstellar reddening for free. For an estimate of 
distances, this is a fundamental advantage of the method. 

In the next step, with effective temperature and gravity measured we can use our synthetic 
spectra to determine metallicity. For this purpose, we concentrate on the objects cooler than 



Urbaneja et al. 



2005a). Three 



17000K, since the S/N is not high enough for the hotter objects (see 
of the cooler objects cannot be used for this purpose, because their effective temperature is not 
constrained by a Balmer jump measurement but by the use of the spectral type already assuming 
solar abundance. In addition, two more objects (Zl 1, CI 6) have spectra too noisy for a metallicity 
fit. Object Z20 shows a metal line spectrum at longer wavelengths, which indicates a spectral type 
somewhat cooler (AO) than the temperature we obtain from the Balmer jump. There is a slight 
chance that this is a composite spectrum, thus, this target is also not used for the determination of 
metallicity (and also not for distance determination, see below). While this reduces the number of 
targets suitable for a metallicity determination, it still leaves us with a sub-sample of 15 objects 
large enough to constrain average metallicity and metallicity gradient of M81, a s we will show 



below . For the measurement of metallicity we apply the technique developed by 



Kudritzki et al 



(2008). For each star we identify spectral windows in the observed spectrum, which are free of 
strong Balmer lines, nebular emission lines or spectral flaws caused by improper correction for 
cosmic ray hits and for which the continuum of the normalized spectrum can be easily matched 
with the one of the synthetic spectra. A pixel-by-pixel comparison of observed and calculated 
normalized fluxes as a function of metallicity then allows for a calculation of ^ 2 ([Z]) in each 
spectral window i and the determination of [Z] ; at which % 2 is minimal. For this comparison, 
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the observed spectra are renormalized for each metallicity so that the synthetic spectrum always 
intersects the qbserv ations at the same value at the edges of the spectral window (see also 



Kudritzki et al 



see again 



2008 ). An average of all [Z], is then used as the measure of metallicity (for details, 



Kudritzki et al. 



2008b . An example is given in Fig[T0]and [T]]for target C20. 



While the analysis method is straightforward and has been tested carefully in previous work, 
two obvious issues, unresolved binarity and blending with fainter sources in the galaxy studied, 
need to be discussed as possible sources of systematic uncertainties. Unresolved binarity can 
affect the analysis in two ways, first, through the contribution of a secondary to the photometric 
fluxes and the spectrum and, second, through the effects of close binary evolution with mass 
transfer or mass loss. In the first case, it is very unlikely that both components have a very similar 
spectral type and luminosity because of the very short lifetime in the supergiant stage. The most 
likely case is a secondary of lower mass still on the main sequence. However, such an object would 
be much fainter by several magnitudes and not affect the spectroscopic analysis or the photometry. 
The second case is more serious, but would affect only the FGLR-distance determination. Binary 
induced mass transfer or mass-loss would change the stellar mass at a gi ven luminosity and c reate 



outl iers from the F GLR-relationship. Such outliers have been found by 



and 



U et al. 



3)- 



Kudritzki et al. 



(2008), 



(2009). They are usually also outliers, when the mass-luminosity relationships of the 



targets are plotted. We will investigate this latter relationship in section 5 (Fig. 14). 

Blending does not appear to b e a problem bec a use o f the enormous optical brightne ss of the 



supergiants as already discussed in 



Kudritzki et al. 



20081) . The study by 



Bresolin et al 



(120051) 



shows that at the distance of NGC 300 at 2 Mpc even ground-based photometry of blue supergiants 
is accurate and not affected by blending. Thus, at 3.5 Mpc for M81 with HST imaging and with 
our careful selection of targets (see section 2) we do not expect blending effects influencing 
the photometry and, therefore, also not the spectroscopy. Of course, in individual cases there is 
always the very small chance of an unresolved coincidence of a target with another bright source. 
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In such cases, the likelyhood that the unresolved blends have the same spectral type is extremely 
small, again because of the short lifetime of blue supergiants. Thus, significant blends should be 
recognized in the spectrum. Target Z20 might be such a case. We also refer the reader to the 
careful m odeling of blending effect s in the HST imaging of Cepheids out to galaxies with 30 Mpc 



distances (Riess et al. 



2009b 



20111) resulting in magnitude corrections of the order of only 0. 1 



mag. Cepheids are 3 to 6 mag fainter than blue supergiants. Thus, since Cepheids are only very 
weakly affected by blending, we do not expect significant effects for supergiants. 

The results of the spectroscopic analysis are summarized in Table 2. Generally, the stellar 
parameters and their unc ertainties are comparable to those obtained in our previous work for 



galaxies less distant (see 



Kudritzki et al. 



2008 



Urbaneja et al. 



2008). We conclude that for this 



type of low resolution quantitative spectroscopy the step from one Mpc (WLM), over 2 Mpc 
(NGC 300) to now 4 Mpc is entirely feasible. In the following, we discuss the results in detail. 



4. Reddening and Extinction 

As described above, one of the advantages of the spectroscopic analysis is that it provides 
information about interstellar reddening. For massive stars imbedded into the dusty disk of a 
star- forming spiral galaxy we expect a wide range of interstellar reddening. Indeed, we find a 
range from E(B-V) = 0.13 to 0.38 mag. Fig [12] shows the distribution of interstellar reddening 
among ou r targets. The a verage value is E(B-V) av = 0.26 mag. The foreground reddening is 
0.08 mag (ISchlegel et al.lll998|) . Our reddening values include both, intrinisc and foreground 
reddening. We stress that our average value of E(B-V) may underestimate the average reddening 
in M81, as our target selection (see Fig. 1) is biased towards lower reddening. 

Fig [12] shows reddening as a function of galactocentric distance. While the scatter is large, it 
is still tempting to fit a regession to the data. We find 
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E(B -V) = (0.415 ±0.025) - (0.0243 ±0.0037) x d/kpc 



(2) 



The lower r eddening beyond 10 



results found by 



Williams et al. 



cpc indicated by this regression is in agreement with the 
(|2009|) . who investigated star formation history and metallicity 
with HST color-magnitude diagrams in the outer fields of M81 and found E(B-V) = 0.14 mag at 
14 kpc galactocentric distance. 

We note that the reddening values found in our study are much larger than the valu e of 



0.03 mag originally assumed in the HST distance scale key project (jFreedman et al. 



1994) for 



Cepheids at inner field s between 3 to 6 kpc galactocentric distance. The final key project study 



(IFreedman et al. 



20011) obtained an average value of E(B-V) = 0.15 mag, still significantly smaller 



than our value, in particular in view of the fact that a difference of 0. 1 mag in reddening results in 
a difference of 0.3 mag in distance modulus if the ratio of total to selective extinction is Ry = 3.1. 



Fig. O (left panel) s 



5. Stellar Properties and Evolution 



rows the location of all t argets in the (log g, log T e yy)-plane compared 



with evolutionary tracks (jMeynet & Maedei 



2003J), which were calculated for solar metallicity 



and which include the effects of rotational mixing and anisotropic mass-loss. The advantage 
of a diagram of this type is that it is independent of any assumption on distance and relies 
completely on the results of the spectroscopic analysis (on the other hand, systematic effects in 
the evolutionary tracks might affect the comparison). The targets form an evolutionary sequence 
crossing from the main sequence towards the red supergiant stage with initial zero age main 
sequence (ZAMS) masses between 15 to 50 M and the majority of objects with ZAMS masses 
about 20 to 25 M . 



A complementary way to discuss stellar evolution and stellar properties is the Hertzprung- 
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Russell diagram (HRD). This requires information about the distance. In section 8 we will use 
the FGLR to determine a distance modulus of jl = 27.7±0.1. With this distance and using the 
spectroscopically determined reddening and extinction and the bolometric corrections provided 
by the model atmospheres for the final parameters of temperature, gravity and metallicity we can 
determine absolute bolometric magnitudes, luminosities and stellar radii. In the calculation of 
stellar radii from luminosities we take into account that the errors in luminosity are dominated 
by the errors in effective temperature and are, thus, correlated (maximum luminosity corresponds 
to maximum temperature and, thus, minimum radius, wheras minimum luminosity at minimum 
temperature yields maximum radius). The results are given in Table 3 and the resulting HRD is 
shown in Fig. [13] (right panel). 



The HRD confirms that the majority of targets is in the ZAMS mass-range of about 20 to 25 
M and is generally consistent with the (log g, log Te/y) -diagram. However, one object (Z15) 
sticks out as very luminous. We recall that the spectroscopic analysis of this object was difficult 
because of extremely strong contamination with nebular H II emission, which might affect the 
determination of gravity in a systematic way which is difficult to assess. In consequence, we have 
not included this object in the FGLR determination of the distance. 

With the stellar radii determined from the luminosities we can use the gravities to estimate 
spectroscopic stellar masses. Those are also given in Table 3. An alternative way to estimate 
masses is to use stellar luminosities and to compare with the luminosities and actual masses at 



the BSG temperatures of evolutionary tracks. Evolutionary masses are al so given in Tab 



Kudritzki et al. 



e 3. 



(2008) 



They are determined from the BSG mass-metallicity relationship given by 
(for Milky Way metallicity and including the effects of rotational mixing). We emphasize that 
both spectroscopic and evolutionary masses are present-day masses and are generally expected 
to be lo wer than the initial Z AMS masses through the effects of mass-loss. Since the early 



work by 



Herrero et al. 



(I1992() it has been found that spectroscopic masses are often significantly 
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smaller than evolutionary masses, although with the development of fully line-blanketed model 
atmospheres and improved NLTE line formation the effect has become much smaller (see 



Kudritzki & Urbaneja 



2009 for a review, and references therein). In Fig. [14] we check our sample 



for this effect by comparing the observed spectroscopic mass-luminosity relationship with the 
prediction of stellar evolution and by directly plotting the ratio of spectroscopic to evolutionary 
mass as a function of luminosity. We find a small effect only at the lower mass end, were 
spectroscopic masses appear to be somewhat smaller than evolutionary masse s. However, we 



concl ude that our sample is not signi ficantly different from the one studied by 



Kudritzki et al 



(2008) in NGC 300 and 



U et al. 



(2009) in M33. 



6. Metallicity, Metallicity Gradient and Chemical Evolution 

Metallicities of 15 targets together with their galactocentric distance are given in Table 2. 
This allows us to discuss stellar metallicity and the metallicity gradient in M81. Fig. [T31 (upper left 
panel) shows a plot of logarithmic metallicity relative to the sun [Z] as a function of galactocentric 
distance (at the distance of 3.47 Mpc - see section 8 - R25 = 11.99 arcmin corresponds to 12.09 
kpc). A metallicity gradient of the young disk populat ion in M81 is cle arly visible. A linear 



regression (using the routine fitexy, Numerical Recipes, 



Press et al. 



1992h yields 



[Z] = (0.286±0.061)-(0.033±0.009)#/fcpc 



(3) 



With respect to the distance independent normalized angular galactocentric distance R/R25 
we obtain 



[Z] = (0.286 ±0.061)- (0.411 ±0. 109) /?//?25 



(4) 
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As is evident from the plot and the regression, young massive stars in the disk of M81 have 
slightly super-solar metallicities at the inner regions and slightly sub-solar metallicity in outer 
parts. The gradient is very shallow, though, compared to the less massive galaxies studied in our 
BSG project. For NGC 30 and M33 metallicity gra dients were det ermined of 0.08 and 0.07 



dexkpc , respectively by 



Kudritzki et al. 



(|2008l) and 



Milky Way, which has a mass comparable to M81, 



U et al 



(2009). On the other hand, for the 



Daflon & Cunhal (120041) in their spectroscopy 



of massive stars obtain a gra dient of -0.031±0.012 dex kpc 1 very similar to our result. (W e 



Rolleston et al 



(2000) forB-stars and 



Luck et al 



note, however, the results by 
Cepheids, who obtained 0.07 dex kpc -1 and 0.055 dex kpc -1 , respectively). 



(2006, 



20 111) for 



Garnett & Shieldsl (|1987|) and lStauffer & Bothunl (|1984|) have analyzed H II region emission 
line spectra of M81 to derive oxygen abundances as a function of ga l actocentric radius. They used 



Pagel et al. 



( 19791) ( Garnett & Shields 



1987 



also 



a strong line method following the calibration by 

used photo-ionization models for an independent check of the abundances obtained). In Fig. 
(right upper panel) we overplot these results with the galactocentric distances corrected to the 
distance used in our work. In the range of 5 kpc to 1 1 kpc there is a large number of objects in 
a similar abundance range as the BSGs with a slight off-set of -0. 1 dex. However, at 5 kpc and 
below there are several objects with very high oxygen abundance. This result might be an artefact 
of the strong-line calibration use d. These inner data point s together with the H II region Muench 1 



19871) lead to an oxygen abundance gradient 



at 16 kpc (carefully discussed in iGarnett & Shieldsl 
of -0.064±0.020 dex kpc -1 with a significantly higher value of [O] = 0.46±0.14 dex at the center, 
where [O] is defined in the same way as [Z], na mely fOl = log (O/Oq) = [O/H] - [O/H] with 



[O/H] = 12 + log (O/H) and [O/H] = 8.69 dex dAUende Prieto et alJl200lh . 



Zaritskv et al. 



(1994) 



The H II regions of M8 1 have also been included in the work by 
who developed a different strong-line calibration method. Their central metallicity is even 
higher, [O] = 0.51±0.11 dex, and the gradient is 0.042±0.015 dex kpc -1 somewhat higher 
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than the result of our BSG work. As has already been shown by iBresolin et al.1 ([2009) and 



Bresolinl (1201 1|) . this calibration leads to metallicities, which are too high when compared with 



H II regio n oxygen abundance s based on electron temperature determinations with auroral lines 



(see also 



Kudritzki et al. 



2008 



for a comparison with BSG metallicities). Our results support this 



conclusion. 



Henry & Howard! (I1995|) used published emission line fluxes of M81 and a series of of 
photoionization models for a study of the oxygen abundance gradient. Their results yield a central 
value of [O] = 0.26 and a gradient of -0.074 dex kpc -1 (Henry, private communication). The 
central value agrees with our BSG work, but the gradient is steeper. 



Stanghellini et al.l (|2010l) have recently studied planetary nebulae (PNe) and H II regions 
in M81 and used the detection of auroral lines to det ermine nebular electron temperatures and 



abundances. Since acc ording to 



a comparison with the 



Bresolin et al 



Stanghe llini et al 



(|2009|) this approach leads to more reliable results, 



(|2010|) H II region oxygen abundances is important. This 
is done in Fig. [15] (lower left panel). At first glance, there seem to be two groups of H II regions, 
one group with abundances comparable to the BSGs and another with abundances 0.4 dex smaller. 
However, for many of the objects t he abundances are too u ncertain with individual errors as 



large as up to 0.6 dex estimated by 



Stanghellini et al. 



(12010 ) and, thus, no clear conclusions are 



possib le with regard to abundance and abun dance gradient from this s ample. 



(I201QQ combine their sample with the one by 



metallicity gradient. However, while the random errors of the 



Garnett & Shield s (1987) to discuss meta 



Stanghellini et al. 



Garnett & Shields 



licity and 



(119871) sample 



are small (0.1 to 0.2 dex), the abundances are affected by the systema tic uncertainties of the 



strong-line method. On the other hand, for the 



Stanghellini et al. 



(2010) abundances the situation 



is opposite, the random errors are large and the systematic errors are strongly reduced. Thus, we 
think the combination of the two samples is subject to uncertainties which are difficult to estimate. 



Contrary to their H II region observations, the PNe analyzed by 



Stanghellini et al 



(12010) 
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have abundances generally more accurate. In Fig. \\5j (lower right panel) they are also compared 
with the BSG metallicities. The average difference in metallicity between the PNe and the BSGs 
is about -0.4 dex and seems to be significant. The metallicity gradient is -0.057±0.007 dex 
kpc -1 and steeper than for the BSGs. This is a very interesting result, since the PNe of this sample 
do not contain type I PNe objects and consist only of type II and III, which means that they are 



signifi 



2010 



cantly older than the BSGs wit h average ages of 3 and 6 Gyrs, respectively (IMaciel et al 



Stanghellini & Haywood 



2010T) . This means that over the last 5 Gyrs the metallicity must 



have increased substantially and the metallicity gradient of the disk has become shallower. 



Williams et al. 



(120091) 



Photometric investigations of the disk of M81 confirm this conclusion, 
in their comprehensive study of star formation and metallicity analyzing HST color-magnitude 
diagrams of an outer disk field at R/R25 = 1.17 find metallicities in the range between [Z] = -0.6 
to -0.3 dex, for the population with ages between 10 Gyrs to 50 Myrs age . They also find solar 



metallicity for the younger population. This result is in ag reement with the 



Tikhonov et al. 



(2005), 



Davidgd (|2009|) . who used the red giant 



who investigated HST CMDs of a different disk field, and 
branch from CFHT MegaCam CMDs over the whole disk of M81 to also estimate a metallicity 
of [Z] =-0.4 dex. While metallicities obtained in this way might suffer from uncertainties in the 
extinction adopted and the systematics of the isochrones used, the picture emerging from the 
combination of our BSG results, the PNe observed and CMDs studied indicates that for a long 
period the metallicity of the M81 disk remained roughly constant and subsolar, but obviously, 
before the birth of the young population of masssive stars, there must have been a phase of 
enrichment. 

This situation is different from the Milky Way. Young massive stars have a metallicity very 



similar to the sun (IPrzvbilla et al. 



sun and to massive stars (Henry et al. 



2008 a|). PNe metallicities are also very clo se to the one of the 



20101 : 



Stanghellini & Haywood 



2010h . The metallicity 



enrichment of the thin disk has been very slow with an estimated increase of metallicity A[Z] 
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0.017 dex Gyr 1 and the metal poor (\Z] -0.58 dex) thick disk may have formed 12 to 13 Gyrs 



ago in a sin gle starburst (IFuh rmann 



Milky Way (IBresolin et al. 



2010 



201 1). We also note that the case of M33 is similar to the 



Urbaneja et al. 



2005b|) . At this point, one can only speculate 



what caused the late en richment of the very y oung population in M81. An interesting thought 



has been formulated by 



Williams et al. 



(12009I). M81 has satellite galaxies such as NGC 3077 



and M82, which are gas and metal ric h (Martin 



M81 (Appl eton et al 



1981 



Heckman et al 



9971) and are involved in tidal interaction with 



1990). Recent inflow from such satellites or the 



tidal interaction induced by them and leading to recent bursts of star formation could then have 
influenced the chemical evolution. 

Chemical evolution models of galaxies also predict changes of the metallicity gradients 
as a function of time, however, many times with qualitatively different results . For instance, 



Chiappini et al. 



(120011) predict gradients to become steeper with time, whereas 



Hou et al. 



(2000) 



predict the opposite . Simulations of disk evolution including the effects of stellar migration by 



Roskar et al. 



(|2008[) also predict a flattening of the gradient through the homogenization of the 



population in the disk as a function of time. 

The comparison of planetary nebulae with a younger stellar generation such as massive stars 
or H II regions offers, in principle, an opportunity to provide observational constraints. In the 



case of the Milky Way 



with time. However, 



Stanghellini & Hav woodl (|2010r) conclude that t he gradient is steepe ning 
J 



Henry et al 



(120101) do 



Maciel & Costal (|2009|) find the opposite, whereas 
not find any hints of evolution at all. Thus, the situation of the temporary evolution of the Milky 
Way abundance gradient re mains controversial. In M 8 1 comparing our BSG results with the PNe 



abundances determined by 



Stanghellini et al. 



(2010) we find a weak indication that the abundance 



gradient became shallower with time. 
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7. Mass - metallicity relationship of galaxies from BSG spectroscopy 



Since the early work by iLequeux et al.l (| 19791) the mass-metallity relationship of star 
forming galaxies has been regarded as an important observational constraint for understanding 
galaxy formation and evolution (see references introduced in the discussion). While these 
pioneering investigations were restricted to a relatively small sample of galaxies, the recent 
spectros copic surveys such as SDSS opened the opportunity to study a large number of such 
objects. iTremonti et al.l (|2004) have analyzed more than 50,000 galaxies observed within SDSS 
and obtained a well defined relationship between oxygen abundance and total stellar mass. 
However, the oxy gen abundance s are a gain based on the use of strong H II region emission 



lines only. While 



Tremonti et al. 



(2004) took special care of this problem and developed their 



own calibrati on of their strong li ne method, the systematic uncertainties are important to be 



investigated. 



Bresolin et al. 



(120091) found that this calibration 



abundances. In a more general approach, 



very 



likely overestimates oxygen 



Kewley & Ellison! (|2008|) demonstrated very clearly 



that the mass -metallicity relationship obtained from the standard strong lines of H II regions 
depends very strongly on the calibration of the strong line method used. Applying ten different 
calibrations, which are frequently used in H II region abundance s tudies , on the same data set of 



emission lines of about 20,000 SDSS galaxies 



Kewley & Ellison 



(120081) obtained the shocking 



result that the mass-metallicity relationship can change from steep to almost flat just dependent on 
the calibration used. Since all the work published with regard to this relationship seems to rely on 
strong line H II region data and given these systematic uncertainties, it seems appropriate to start 
an investigation based on stellar spectroscopy only. With the results obtained here and compiling 
the metallicities of the BSG quantitative spectroscopy work for other galaxies published so far we 
have made a first attempt. 

The compilation of galaxy masses and metallicities is given in Table 4. For the spiral galaxies 
with a clear metallicity gradient (NGC300, M33, MW, M31, M81) metallicity values were taken 
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at galactocentric distances of two disk scale lengths. For the irregular Local Group galaxies 
average values were used. The data are plotted in Fig. [16] (left panel). A very clear correlation of 
metallicity with stellar mass is obtained. 

While the weakness of our approach at this stage is the small size of our sample, it is 
tempting to compare with the SDSS H II region based results discusse d. For this purpose we 
have overplotted the average mass -metallicity relationships obtained by iKewley & Ellison! (12008) 



for the ten different calibrations used in their work. It seems that a few of these calibrations 

i ii 1 i 



Tremonti et al.l 120 04: Zaritskv etal 



whereas others (IPettini & Pagel 



1994 ) lead to a much steeper relationship than our work, 
|2004|) are in much better agreement. We note that our sample is 
probing a larger galaxy mass range than the SD SS studie s , goin g from low-mass dwarf irregulars 



to giant spirals. As pointed out in the study by 



Lee et al. 



(|2006|) this is important for 



the scenarios for galaxy formation and evolution. (We realize that in 



Lee et al. 



constraining 



2006 the stellar 



masses of some of the d warf irregu 



than the masses given by 



Woo et al 



ars o verlapping with our sample are significantly smaller 



2008 



which we use for Fig. [161 This will require further 
investigation). In future work we plan to enlarge the sample of galaxies with quantitative studies 
of BSGs to make this comparison more significant. 



8. Distance 



The FGLR is a tight correlation between the flux- weighted gravity (gp = g/^ 4 eff, Teff in 
units of 10 4 K) and the absolute bolome tric magnitude M& / of BA supergiants. As described 



in detail in iKudritzki et al 



(2003, 



2008) the physical background for this relationship is the fact 



that massive stars evolve at constant luminosity and mass accross the HRD from the hot main 
sequence to the red supergiant stage. During this evolution, gf remains constant, because of 
the constant luminosity and mass. On the other hand, stellar luminosity is a strong function of 
stellar mass (see Fig. [14] as an example) and, therefore, also a strong function of flux- weighted 
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gravity , which establishes the FG LR. (For all det ails, we refer the reader to the two papers just 



cited). 



Urbaneja et al 



(120081) and 



U et al 



(2009) were the first to use the FGLR for distance 



determination of the metal poor dwarf galaxy WLM and M33, respectively. Here, we follow the 
same procedure as was detailed in these papers. 



The FGLR has the form 



M bol = a(\ogg F -l.S) + b 



(5) 



with the recent calibration provided by lKudritzki et al. 



(2008), a = 3.41 and b = -8.02. 



For each of our targets the spectroscopic analysis yields de-reddened apparent bolometric 
magnitude m^i and flux-weighted gravity, which are given in Table 2. These data are plotted in 
Fig-ELZl Very obviously, there is a clear relationship between flux-weigthed gravity and apparent 
bolometric magnitude. We can use these data to fit a regression of the form 



m bo i = a(\ogg F - 1.5) + Z?M8i • 



(6) 



The fit resul t is also shown in Fi g. [171 Since our targets span only a limited range in gp 



compared to the iKudritzki et al. 



(2008) calibration sample, we adopt the slope value provided 
by this calibration and fit only the intercept &M81- The difference between b and £>M8i yields 
the distance modulus, which we d etermine to be /i = 27.71 ±0.08 mag (the error is calculated 



similarly as in 



Urbaneja et al 



2008). 



The 



Kudritzki et al. 



(2008) calibration of the FGLR is based on data from eight galaxies 
with distances mostly determined from using Cepheids. Recently, we have started the study of 
a large sample of BA supergiants in the LMC using high resolution, high S/N spectra with the 
goal to provide a new calibration of the FGLR based on the LMC only. This work is almost 
completed and will be published soon (Urbaneja et al. 2011, to be submitted to ApJ). With an 



adopted distance modulus to the LMC of m-M = 18.50 mag we obtain the calibration values &lmc 
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= 4.53 and \>lmc = -7.88. While this is a significantly steeper FGLR at the low luminosity /high gf 
end, this change in calibration does barely affect our distance determination, because most of our 
targets are at lower g/7 /higher luminosity. A regression fit with these (still preliminary) calibration 
values yields a distance modulus of /i = 27.68 ±0.09 mag. We, thus, adopt a distance modulus of 
/I = 27.7 ±0.1 mag. 

We compare this value with previ ous distance determinations b ased on Cepheids. In addition 



to the HST Key Project wor k on M8 1 (IFreedman et al. 



by 



McCommas et al 



the Wesenheit method (Madore 



1994 



2001) there are two recent studies 



(200 91) and bvlGerke et al] (1201 Cepheid distance studies typically apply 



1982|) with a combination of V and I band magnitudes which 
is assumed to be reddening free and then compare with th e corre sponding period luminosity 



relationship of LMC Cepheids. Following 



Kennicutt et al. 



(1998), distances are corrected for 



the difference in abundance between the target Cepheids and those in the LMC. This so-called 
"metallicity correction" has the form Aji = y([0/H] - 8.5) where [O/H] = 12 + log (O/H) is 
the logarithmic oxygen abundance of the young stellar population in the target galaxy at the 
galactocentric distance of the observed Cephei d field relative to hydrogen. 7 is a fit parameter and 



has been determined by 



Kennicutt et al. 



(11998) from the fact that Cepheids in inner fields of the 
spiral galaxy M101 are brighter and yield a shorter apparent distance modulus than those in outer 
fields. Attributing this difference to a metallicity dependence of the period luminosity relationship 



and adopti ng stellar meta l 



studies by 



Zarit skv et al 



icities and metallicity gradien ts from the oxygen H II region strong line 



dl994h . 



Kennicutt et al 



(119981) obtained 7 = -0.29 mag dex -1 . [O/H] 



= 8.5 dex in this metallicity correction is the adopted value of this abundance for the LMC. It 
refers to the "old" oxygen abundance scale where where [0/H] Q =8.9 dex. (We will show below 



that this va^ 



Macrietal 



ue is to o high independent of the actual value of the oxygen abundance for the sun). 



(|2006l) found a similar value of 7 for the maser galaxy NGC 4258 again from the 



different distance moduli obtained from inner and outer field Cepheids. 
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McCommas et al.1 (|2009h in their Cepheid distance investigation of M8 1 use HST light curves 
of 1 1 fundamental and two first overtone short period Cepheids in the outer disk of M8 1 at R 
= 1.23 R25 (~ = 13.5 kpc) and obtain a distance modulus of M81 relative to the LMC of A/i = 
9.34 ± 0.05 mag. Checking the consistency with the 25 long period Cepheids i n two inner HST 



WFPC fields observed by the Key Project located at R = 0.36 R25 (~ = 4.3 kpc) 



McCommas et al 



(|2009|) use the same Wesenheit forma l ism a nd obta in a distance mo dulus 0.23 mag shorter. 



Following the work by 



Kennicutt et al. 



Macri et al 



(12006 ) they also apply a metallicity 



(119981) and 

correction with 7 = -0.29 mag dex -1 . This correction introduces a small increase of the distance 
to A/i = 9.37 ±0.05 mag and reduces the difference in distance mod ulus between outer a nd inner 
field Cepheids to 0.09 mag. It is based on the metallicity study by IZaritsky et al.l (| 19941) who 
obtained [O/H] = 9.196 - 0.49 R/R25 for the oxygen abundance as a function of galactocentric 
distance as a result of their strong-line analysis of H II region emission lines. Explaining the full 
diffe rence in distance mo dulus between inner and outer field Cepheids in terms of metallicity with 



the 



Zari tskv et al 



(|1994|) metallicity gradient requires 7 = -0.55 mag dex 



-1 



Gerke et al.l (|201lf) investigate 107 long period Cepheids observed with the LBT in a 
galactocentric range of 0.29 < R/R25 < 0.88 and with ground-based B, V, I photometry. Without 
applying a metallicity correction they obtain A/i = 9. 19 ±0.05 mag. They also realize a trend 
in Cepheid distance modulus as a function of galactocentric distance and obtain a metallicity 
correction, which leads to 7 = -0.56 ±0.36 mag dex" 1 and a distance modulus of A/i = 9.39 ±0.14 



mag. This agrees with with 



McCommas et al 



(|2009|) and also with the original value of the Key 



Project of A^i = 9.30 ±0.15 mag 

Our FGLR distance to M81 is based on a LMC distance modulus of 18.5 mag and, thus, 
a difference of A/i = 9.2 ±0.1 mag. This is 0.10 to 0.19 mag or 5 to 8% shorter than the ones 
obtained with the Cepheid work. However, we note that there is good agreement with the inner 
field long period Cepheids, when no metallicity corrections are applied. In the following we 
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discuss some aspects of this metallicity correction. 



With the solar o xygen abundance [O/H] = 8.69 dex (lAllende Prieto et al. 



Zaritsky et al. 



200 lh the 

(11994)) logarithmic oxygen abundances relative to the sun are [O] = [O/H] - [O/H] 



= 0.506 - 0.49 R/R25- If oxygen is taken as proxy for metallicity, this is a significantly higher 
metallicity than found in our BSG spectroscopy in equation (4), while our gradient is shallower. 
Applying our metallicity gradient to correct for distance modulus difference between the inner 
and outer field Cepheids in M81 would require an even more negative value of y, namely 7 = -0.65 
mag dex -1 . Moreover, the LMC oxygen abundance [0/H]lmc = 8.50 dex or [0]lmc = -0.19 dex 
adopted i n these corre c tions is too large compared with the LMC oxygen abundance of B-stars 



found by 



Hunter et al. 



(|2007|) (\0 /H]i,mc = 8.33 dex or \0) j,mc = -0.36 dex ) ? the i ron abundances 



of LMC Cepheids determined by 



Romaniello et al. 



(2008D and 



Luck 



et al. 



1998J) (f 



Bresolinl d201lh ([0/H] LMC 



-0.33 dex), and the LMC H II region oxygen abundances obtained by 
= 8.36 dex or \0\lmc = -0.33 dex). This means that with our BSG metallicity values in M81 the 
Cepheids in the outer field have a metallicity 0.1 1 dex higher than the LMC. If one would apply 
the metallicity correction with 7 = -0.65 mag dex -1 accordingly, this would enlarge the distance 
modulus by another 0.07 mag. 

However, with such a large negative value of 7 it is important to note that this empirical 
correction for the metallicity dependence of the period-luminosity relationship, which claims that 
Cepheids become brighter with increasing metallicity, is in striking disagreement with pulsation 



theory, which predicts exactl y the opposite, name 



with increasing m etallicity (jFiorentino et al 



Bono et al. 



2002; 



y that the Cepheid brightness decrea ses 



Marconi et al. 



2005 



Fiorentino et al. 



2007; 



2008J). It also disagrees with the recent high S/N, 



ligh spectral reso l ution q uantitative 



Romaniello et al. 



(2008), which 



spectroscopy in the Milky Way and the LMC carried out by 
confirms the prediction by pulsation theory. According to this work, the value of 7 should be 
positive and not negative. In other words, as careful spectroscopic metallicity studies compared 
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with observed differences of distance moduli between inner and outer field Cepheids push 7 to 
increasingly negative values, an explanation of that distance modulus differences in terms of 
metallicity seems unlikely. It must be something else and it is an additional systematic effect not 
understood. 



We also note that lU et al.1 ((2009) have demonstrated from their quantitative spectroscopy 
of blue supergiants in M3 3 that the diffe r ence o f distance moduli between inner field and outer 



field Ce pheids found 



worse, iBresolin et al. 



by 



Scowcroft et al 



(|2009() would require a y-value of -0.55 mag/dex. Even 
(|2010|) re-determined H II region abundances in M33 using auroral lines and 
applying thei r abundance gr adient to the Cepheid fields in M33 yields 7 = - 1.2 mag dex -1 (see 



discussion in 



Bresolin 



2011 ) 



Another galaxy where the comparison of Cepheids in the inner and outer fields leads to 
a significantly different distance modulus is the maser galaxy NGC 4258. This galaxy is of 
particula r importance, since it has b een used as the new anchor point for the extragalactic distance 



scale by 



Riess et al. 



(20093^, 



20111 ) because of its accur ately known distance fro m the Keplerian 



motion of water m asers orbiting the central black hole ([Humphreys et al. 



Macri et al. 



2008). However, 



(|2006h . who carried out the HST obervations of Cepheids in NGC 4258 again found 



the distance modulus of the inner field Cepheids to be sho rter than in the oute r fields and based on 



the H II region strong line method o xygen abundanc es by 



of -0.29 mag dex . Most recently, 



Zaritsky et al 



1994) derived a 7- value 



Bresolinl (|201ll) re-determined the H II region metallicities 



in this galaxy including the observation of auroral lines in a few cases. This led to a downward 
substantial revision of the metallicity, which seems to be close to the LMC and not strongly 



Bresolin (2011) show 



super-solar, and a very shallow abundance gradient. Based on these results, 
that 7 = - 0.69 mag dex -1 would be needed to explain the distance modulus difference between 
inner and outer fields, again a value much too negative, when compared with pulsation theory and 
observational work on Milky Way and LMC Cepheids. While the improved H II region work on 
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this important galaxy still awaits an independent confirmation through a study of BSGs, it is an 



this p oint. 



Majaess et al. 



(|201ll) discuss the large metallicity co rrections suggested b y 



(1201 Y) and by the recent HST/ACS Cepheid study of M 101 by 



Shappee & Stanek 



v Ger 


ce et al. 


(2011 


) and 



demonst rate that such corre ctions lead to very improbable distances of the LMC and SMC. The 



work 



by 



Storm et al 



Majaess et al. 



(|201l|) indicates that a lower limit for 7 is -0.2 mag dex -1 . 
(|2011|) argue that crowding is very likely responsible for the distance modulus differences obtained 
between inner and outer field Cepheids and not metallicity. We think that a careful spectroscopic 
investigation of galactic metallicities and their gradients and distance determinations using the 
FGLR as an independent method will help to clarify the situation. 

Independent of the Cepheid work there have been numerous studies of HST color-magnitude 



diagrams of M8 1 to determine a distance from t he tip of the red giants branc 



moduli found wer e 28.03 mag (|Sakai et al 



(Rizzi et al. 



20071) . 27.72 to 27.78 mag (IDalcanton et al 



1. The distance 



200 41), 27.93 



2009 



mag (ITikhonov et al 



20051) . 27.70 mag 



different fields in the halo and 



the outer di sk), 27.81 mag (E xtragalactic Distance Database catalogue, 



Tullv et al. 



27.86 mag (IDurrell et al. 



2009) and 



I2010I) . The more recent work since 2007 has converged on an improved 
methodology and seems to agree, within the uncertainties, with the distance modulus found in our 
study. 



9. Conclusions and Future Work 

In this paper we have demonstrated that the quantitative spectroscopy of BSGs is a promising 
tool to constrain the chemical evolution of galaxies and to determine their distances, which can be 
applied to galaxies clearly beyond the Local Group. Using the relationship between flux-weighted 
gravity and luminosity we were able to determine a new distance to M81, which compares well 
with TRGB distances. While there is also agreement with HST Cepheid distances within the 
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error margins, our results with regard to metallicity and metallicity gradient confirmed previous 
studies that the systematic differences between distance moduli obtained from inner and outer 
field Cepheids (found in M33, M81, M101, NGC 4258) are very likely not caused by a metallicity 
dependence of the period-luminosity relationship of Cepheids. There must be another reason for 
these systematic differences. 

An independent check of distances obtained with either the TRGB or Cepheids is important 
for future work. We note that besides the importance for characterizing the physics of galaxies in 
the Local Volume accurate distances and a careful discussion of the systematics of stellar distance 
determination methods are cr ucial for constrain ing the dark energy equation-of-state parameter w 



= p/(pc 2 ). As is well known (IMacri et al. 



2006), the determination of cosmological parameters 
from the cosmic microwave background is affected by degeneracies in parameter space and 



cannot provide strong constraints on the value of Hq (ISpergel 



2006 



Teg mark et al. 



20041) . Only 



if additional assumptions are made, for instance that the universe is flat, Ho can be predicted with 
high precision (i.e. 2%) from the observations of the cosmic microwave background, baryonic 
acoustic oscillations and type I high redshift s upernovae. If these assumptions are rela xed, 



then much larger uncertainties are introduced (|Spergel et al.l 



2007; 



Komatsu et al. 



2009). The 



uncertainty of the determination of w is related to the uncertainty of Ho through Aw/w pa 2AHo/Ho. 
Thus, an independent determination of Ho with an accuracy of 5% will allow the uncertainty 
of w to be reduced to 0.1. While extremely p ro mising steps towards this goal have been made 



by 



Macri et al 



(2006() and 



Riess et al 



(2009a,b 



20111) using the maser galaxy NGC 4258 as a 



new anchor point and HST IR Cepheid photometry of recent SNIa galaxies out to 30 Mpc, it 
is clear that the complexity of this approach requires additional and independent tests. Crucial 
contributions which can be made using BSGs besides independent distance determinations are to 
investigate the role of metallicity and interstellar extinction. 



We have also shown that the determination of metallicities for individual supergiant stars 
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beyond the Local Group is possible. In this way, we can determine galaxy metallicities and 
metallicity gradients avoiding the systematic uncertainties of H II region strong line methods. 
This can be used as an independent way to directly measure the mass -metallicity relationship 
of galaxies and to correlate metallicity gradients with galactic properties such as mass, angular 
momentum and morphological type. But it can also be used to find out about systematic 
uncertainties of H II region strong line method calibrations and to identify the more reliable 
ones or to develop a new one tested with BSG metallicities. Moreover, in combination with 
metallicity information of an older population of stars obtained through the analysis of CMDs or 
the spectroscopy of PNe the chemical evolution history of galaxies can be investigated. In the case 
of the disk of M81 we have found an indication of a late enrichment of heavy elements, which 
is significantly different from the Milky Way. We have also provided the first mass-metallicity 
relationship for star forming galaxies solely based on stellar spectroscopy. 
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Fig. 1. — Selection of M81 BSG targets. Left: Color magnitude diagram (photometry from 
Dalcanton et al.ll2009Q with selection box (blue dashed) and selected targets (red). Right: Location 
of selected targets within M8 1 . 
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Fig. 2. — Analysis of object Z10. Left: Fit of the observed Balmer jump. The final model with the 
parameters given in Table 2 (thick solid curve) fits the Balmer discontinuity well. Two models with 
T e ff higher/lower by 500K are also shown (dashed) to demonstrate the temperature sensitivity of 
the D s fit. Right: Fit of two Balmer lines with the final model. Two models with log g higher/lower 
by 0.1 dex are shown (dashed) to demonstrate the gravity sensitivity of the Balmer line fits. Note 
that the strong spectral line at the left edge of the panel for Hio is Hn, which is not used for the 
fits, because it is at the edge of the normalized spectrum, where continuum rectification becomes 
difficult. 




Fig. 3. — Same as Fig. 2 but for object C20. 
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Fig. 4. — Fit diagrams for the fit of the Balmer jump (steeper curves) and the Balmer lines. Left: 
the (log g, T e ff) diagram, right: (log g^, T eff ). Teff is given in 10 4 K. The dashed curves indicate 
maximum errors of the fits. For discussion, see text. 




Fig. 5. — Balmer jump fit for 12 objects in field Z. Logarithm of flux is plotted vs. wavelength in 
A. The bar in each panel indicates 0.05 dex changes in flux level. 
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Fig. 7. — Balmer jump fit for 6 objects in field C. Logarithm of flux is plotted vs. wavelength in 
A. The bar in each panel indicates 0.05 dex changes in flux level. 
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Fig. 8. — Balmer line fit for 6 objects in field C. Normalized flux is plotted vs. wavelength dis- 
placement from the line center in A. 
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Fig. 9. — Si III line fits (left), Si IV (solid bars) and Si II (dotted bars) line fits (middle), and 
hydrogen H10 line fits (right) for the three early B supergiants of our sample. Note that Si IV 
41 16A is blended by He I. Normalized flux is plotted vs. wavelength in A. 
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Fig. 10. — Metallicity fits in 6 spectral windows of object C20. The synthetic spectra are plotted 
in bold and the metallicity is indicated at the left bottom of each plot. Normalized flux is plotted 
vs. wavelength in A. Fe lines are indicated by solid bars, Cr: dashed, Ti: dotted, Si: dash-dotted, 
Mg: dash-triple dotted. Only 6 of the 14 available metallicities are shown. 
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Fig. 11. — # 2 ([Z]) for each spectral window of object C20 as a function of metallicity [Z]. The 
curve for each window has a well defined minimum abcissa [Z] ; . The average of all [Z], is adopted 
as the stellar metallicity value. 



-50- 




0.00 0.10 



0.20 0.30 
E(B-V) 



0.40 




Fig. 12. — Interstellar reddening in M81. Left: Histogram of the E(B-V) distribution. Right: E(B- 
V) as a function of galactocentric distance with the regression curve (dashed) discussed in the text.. 




Fig. 13. — Stellar parameters of the observed sample of M81 supergiants compared with 
evolutionary tracks for the Milky Way metallicity including the effects of rotational mixing 



( Mevnet & Maeder 



2003b . Left: (log g, log T e ff) - diagram. Right: Hertzsprung-Russel diagram. 



The zero-age main sequence masses are (in increasing luminosity /decreasing gravity) 12, 15, 20, 
25, 40 solar masses, respectively. 




Fig. 14. — Left: Observed mass-luminosity relationship compared with stellar evolution theory 
using the tracks from Fig. [13] at an effective temperature of 10 4 K. Right: Logarithmic ratio of 
spectroscopic to evolutionary masses as a function of luminosity. 
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Fig. 15. — Upper left: Metallicity of BSGs in M81 as function of galactocentric distance in kpc. 
The dashed curve is the regression discussed in the text. Uncertainties are given in Table 2 and 
not plotted. Upper right: S ame as upper left with the the oxygen abundances of H II regions from 



the strong line studies by 



Garnett & Shields 



Stauffer & Bothunl (11984) overplotted. 



(11987b and 

Random uncertainties of the H II region data are between 0,1 to 0.2 dex, systematic uncertainties 
ar e discussed in the text. Lo wer left: Same as upper left but with the H II region oxygen abundances 



by 



Stanghellini et al 



of PNe obtained by 



(20 101) overplotted. Low er right: Same as upper left, but oxygen abundances 



Stanghellini et al. 



( 20101) overploted with error bars. The dashed line in all 



four panels is the BSG regression obtained in this work. For a detailed discussion, see text. 
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Fig. 16. — Left: Observed mass-metallicity relationship of galaxies obtained from spectroscopic 
studies of blue supergiants. The red square is the M81 result from this paper. Right: Same as 
left, but now with the average relationships obtained by iKewley & Ellison! (120081) for the ten dif- 
ferent H II region strong lin e calibrations u s ed in their study of 20,000 SPSS galaxie s. The ten 



cali brations are: |T1 (solid ) 



tet) 



Kobu 



McGaugh 



nicky 



Tremonti et al 



(2004), 



& Kewlevl (120041) . [4] ( dash-dotted' 



21 (dashed) 



Kew 



19911) . [6] (dash-triple-dotted) iDenicolo et al. 



(using TO 



III1/H ^ and rNHl/H,y ), T81 (dashed) 



Pilyuginl (120011) . [10] (dotted) 



Zaritskv et al 
ij 



1994), [3] (dot- 



ev & Dopital (120021). T51 (long-dashed) 

4 



(2002), [7] (solid lPettini & Pagell 



Pilyugin & T huan (2005). 



Pett 
iaJ ( 



(2004 



ni& Pagell (12004) (using [Nll]/H a ), [9] (dotted) 
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Fig. 17. — The observed FGLR in M81. Solid squares are targets used for the distance determi- 
nation fit. Targets plotted as open squares were not included in the fit for reasons ex plained in the 



Kudritzki et al. 



(2008). The dashed- 



text. The dashed line corresponds to the FGLR calibration by 
dotted line is the new (still preliminary) LMC calibration (Urbaneja et al, 201 1, to be submitted to 
ApJ) discussed in the text. Both calibrations yield a very similar distance modulus. 
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Table 1. M81 - Spectroscopic targets 



No. 


name 


«2000 


52000 


R/R 25 


sp.t. 


my 


B-V 










h min sec 


o / // 






mag 


mag 


dex 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 






Zl 


9 55 30.580 


69 12 22.716 


0.83 


B7 


20.946 


0.052 


0.126 


1 


Z2 


9 55 34.965 


69 11 58.524 


0.81 


B4 


21.493 


0.011 


0.084 


2 


Z3 


9 55 25.022 


69 11 16.908 


0.69 


B7 


20.847 


0.093 


0.122 


3 


Z4 


9 55 12.686 


69 10 50.880 


0.61 


B9 


20.305 


0.198 




4 


Z5 


9 55 26.224 


69 10 16.896 


0.60 


B3 


21.452 


0.059 


0.056 


5 


Z6 


9 55 18.204 


69 09 51.984 


0.53 


B7 


21.206 


0.185 


0.156 


6 


Z7 


9 55 30.948 


69 09 33.696 


0.55 


A3 


21.454 


0.244 


0.539 


7 


Z9 


9 55 10.488 


69 08 27.132 


0.41 


A4 


21.296 


0.358 




8 


Z10 


9 55 34.353 


69 08 39.552 


0.48 


B9 


21.236 


0.185 


0.230 


9 


Zll 


9 55 21.518 


69 08 15.864 


0.39 


Bl 


21.327 


0.048 




10 


Z12 


9 55 35.100 


69 08 16.908 


0.45 


B4 


21.020 


0.185 


0.049 


11 


Z13 


9 55 29.344 


69 07 48.432 


0.37 


B1.5 


21.113 


0.113 




12 


Z14 


9 55 31.761 


69 07 39.036 


0.36 


B7 


21.371 


0.207 


0.127 


13 


Z15 


9 55 34.783 


69 07 31.440 


0.37 


B0.5 


20.495 


0.136 




14 


Z16 


9 55 43.579 


69 07 18.768 


0.42 


B4 


19.979 


0.135 


0.002 


15 


Z17 


9 55 39.237 


69 06 35.172 


0.69 


A4 


21.165 


0.372 


0.586 


16 


Z18 


9 55 43.442 


69 06 18.972 


0.33 


B9 


20.402 


0.211 


0.171 


17 


Z20 


9 55 46.972 


69 05 48.516 


0.32 


Al 


20.330 


0.305 


0.110 


18 


C6 


9 54 35.976 


69 05 00.168 


0.74 


Al 


20.784 


0.306 


0.288 


19 


C9 


9 54 51.542 


69 05 33.288 


0.51 


B6 


21.217 


0.071 


0.129 


20 


Cll 


9 54 49.214 


69 06 17.640 


0.53 


B9 


20.412 


0.296 


0.178 


21 


C13 


9 54 36.451 


69 07 15.708 


0.69 


B2 


21.152 


0.038 


0.015 


22 


C14 


9 54 35.196 


69 07 43.752 


0.70 


B3 


21.414 


0.039 


0.055 


23 


C16 


9 54 54.530 


69 08 14.892 


0.51 


B9 


21.245 


0.337 


0.186 


24 


C20 


9 54 51.079 


69 09 43.992 


0.60 


B9 


20.411 


0.259 


0.249 


25 


C21 


9 55 18.777 


69 09 52.668 


0.53 


B8 


20.395 


0.105 





b 



b 



b 



Table 1 — Continued 



No. name a 2 ooo ^2000 R/R25 sp.t. my B-V D B 
h min sec ' " mag mag dex 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) 



"Galactocentric distance, in units of R25 = 11.99 arcmin ~ 12.09 kpc (distance modulus 
27.70 mag). A position angle PA = 157°, an inclination i =57°and central coordinates 
«2000 = 9h55min33.2sec, §2000 = 69°3'55"were assumed (Hyperleda data base, Paturel et 
al., 2003) 

b no near UV spectral coverage; no Balmer jump measured 
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Table 2. Stellar Parameters 



No. 


name 


T eff 


logg 


log gf 


[Z] 


E(B-V) 


BC 










K 


cgs 


cgs 


dex 


mag 


mag 


mag 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 






Zl 


1 2500 1020 


1 75° 17 
1 -' J 0.20 


1 36° 05 
i " 5o 0.06 


1 


Z2 


15000™ 


2.20°;" 


1 50° 05 
1 " 5U 0.05 


2 


Z3 


12500?«° 


l - /J 0.21 


1 34O.O5 
1 -'^0.06 


3 


Z4 


10000!°° 


1 45 a 17 


1 45O.O8 


4 


Z5 


16000?* 


2 15 a 14 
Z - 1J 0.17 


1 o->0.05 
1 - J -'0.05 


5 


Z6 


12500"° 


1 95 a 15 




6 


Z7 


8500!«° 


1 40°- 17 


1.688;!? 


7 


Z9 


8300|°° 


1 20°- 19 


1 52° 15 
1 - JZ 0.20 


8 


Z10 


11000^° 


1 75°' 13 
i. u Q 14 


1 -- )8 0.07 


9 


Zll 


22000'°°° 
zzuuu 1000 


2 62° 18 
Z ' DZ 0.18 


1 25°' 10 


10 


Z12 


15000jl°° 


1 95 a 17 


1 25°'° 5 


11 


Z13 


2I000 1000 
zluuu 1000 


? S4O.I8 
Z - J T).18 


1 25O.IO 


12 


Z14 


13000™ 


1 90°- 15 
lyu 0.17 


1 440.05 


13 


Z15 


25000 1000 


2.64 


1.05 


14 


Z16 


, <ooo 10 °° 

13uuu 1000 


1.80°;i« 


1 10°° 5 

1 - lu 0.05 


15 


Z17 


8300'i° 


, O.-0.14 
1 -' J 0.18 


1 67°' 12 


16 


Z18 


11500™ 


1 65°- 17 


1 41 0.06 


17 


Z20 


12000!'«° 


1 55°- 23 
1 - J - , 0.27 


1 2^0.05 

'•"o.os 


18 


C6 


9250«° 


1.20°;?1 


1 34°- 11 
'•-"0.09 


19 


C9 


13000™ 


1.908:!? 


1 44O.O5 


20 


Cll 


11000"° 


1 55° 17 


1 ISO.06 


21 


C13 


17000™ 


2 270.13 
z - z '0.15 


1 35°° 5 
1 ' JJ 0.05 


22 


C14 


17000™° 


2 07O.I2 
Z - J/ 0.13 


1 45O.O5 
1 ' 4J 0.05 


23 


C16 


11000™ 


1 55°' 17 
l - JJ .19 


Li ' i o.m 


24 


C20 


10500™, 


1 60°- 15 


1 coo- 07 
1 - JZ 0.08 


25 


C21 


12500!°° 


1 75° 12 


1 3fi°° 5 



-0.07±0.15 


0.13 


-0.75 


19.78±0.18 




-0.01±0.2() 


0.15 


-1.15 


19.88±0.16 




-0.09±().1() 


0.17 


-0.75 


19.56±0.19 






0.22 


-0.28 


19.36±0.13 


a 


0.03±0.20 


0.21 


-1.29 


19.54±0.17 


b 


0.08±0.10 


0.28 


-0.73 


19.61±0.15 




-0.03±().2() 


0.23 


0.02 


20.76±0.10 






0.31 


0.05 


20.35±0.10 


a 


0.11±0.10 


0.25 


-0.45 


20.03±0.13 






0.24 


-2.12 


18.47±0.15 


c 


0.07±0.10 


0.30 


-1.17 


18.92±0.20 






0.29 


-2.01 


18.19±0.16 


c 


0.21±0.15 


0.31 


-0.83 


19.59±0.17 






0.32 


-2.45 


17.04±0.17 


c 


0.09±0.10 


0.27 


-1.29 


18.08±0.17 


b 


0.14±0.15 


0.34 


0.05 


20.13±0.10 




0.1 5 ±0.10 


0.27 


-0.56 


18.99±0.17 






0.36 


-0.69 


18.51±0.23 




0.00±0.10 


0.28 


-0.16 


19.75±0.23 






0.17 


-0.83 


19.86±0.17 




0.04±0.15 


0.34 


-0.47 


18.88±0.15 






0.19 


-1.49 


19.09±0.16 


b 




0.19 


-1.49 


19.33±0.15 


b 




0.38 


-0.47 


19.58±0.15 




0.08±0.10 


0.30 


-0.36 


19.12±0.15 






0.19 


-0.75 


19.07±0.13 


a 
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Table 2 — Continued 



No. name T e g log g log g F 
K cgs cgs 

(1) (2) (3) (4) (5) 



[Z] E(B-V) BC m bo , 
dex mag mag mag 

(6) (7) (8) (9) (10) 



"no Dg, T e g from spectral type 
b 7^from Dg and Si II, Si III, Si IV 
c 7^from Sill, Si III, Si IV 
d extreme HII contamination of Balmer lines 
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Table 3. Absolute magnitudes, luminosities, radii and masses 



No. 


name 


M v 




log L/L a 


R 




Meval 






mag 


mag 


dex 


R s 


M e 


M e 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 






Zl 


-7.16 


-7.92 


5.07±0.07 


73.1±6.2 


10.9 


17.0 


1 


Z2 


-6.67 


-7.82 


5.03±0.06 


48.5±3.5 


13.5 


16.4 


2 


Z3 


-7.38 


-8.14 


5.16±0.08 


80.9±7.1 


12.6 


18.2 


3 


Z4 


-8.08 


-8.34 


5.24±0.05 


138.6±8.4 


19.7 


18.2 


4 


Z5 


-6.90 


-8.16 


5.16±0.07 


49.8±3.9 


12.7 


18.3 


5 


Z6 


-7.36 


-8.09 


5.14±0.06 


79.1±5.6 


20.2 


17.9 


6 


Z7 


-6.96 


-6.94 


4.68±0.04 


100.7±4.9 


9.2 


12.8 


7 


Z9 


-7.37 


-7.35 


4.84±0.04 


127.6±6.1 


9.4 


14.3 


8 


Z10 


-7.24 


-7.67 


4.97±0.05 


84.2±4.9 


14.5 


15.7 


9 


Zll 


-7.12 


-9.23 


5.59±0.06 


43.2±3.0 


28.2 


27.0 


10 


Zll 


-7.61 


-8.78 


5.41±0.08 


75.5±6.9 


18.4 


22.4 


11 


Z13 


-7.49 


-9.51 


5.70±0.06 


53.9±3.9 


36.5 


30.1 


12 


Z14 


-7.29 


-8.11 


5.14±0.07 


73.8±3.9 


15.7 


18.0 


13 


Z15 


-8.20 


-10.66 


6.16±0.04 


64.6±3.3 


66.1 


52.6 


14 


Z16 


-8.56 


-9.62 


5.75±0.07 


111.1±8.6 


28.3 


30.7 


15 


Z17 


-7.59 


-7.57 


4.93±0.04 


141.2±6.6 


16.2 


15.3 


16 


Z18 


-8.14 


-8.71 


5.38±0.07 


124.3±9.6 


25.1 


21.9 


17 


Z20 


-8.49 


-9.19 


5.58±0.09 


142.4±15.3 


26.1 


26.0 


18 


C6 


-7.78 


-7.95 


5.08±0.05 


135.4±8.4 


10.5 


17.1 


19 


C9 


-7.01 


-7.84 


5.04±0.07 


65.2±5.0 


12.2 


16.5 


20 


Cll 


-8.34 


-8.82 


5.43±0.06 


142.9±10.1 


26.3 


22.7 


21 


C13 


-7.14 


-8.61 


5.34±0.06 


54.3±4.0 


19.9 


21.2 


22 


C14 


-6.88 


-8.37 


5.25±0.06 


48.6±3.3 


20.1 


19.5 


23 


C16 


-7.63 


-8.12 


5.15±0.06 


103.5±7.3 


13.8 


18.0 


24 


C20 


-8.22 


-8.58 


5.33±0.05 


140.4±8.4 


28.5 


21.0 


25 


C21 


-7.89 


-8.63 


5.35±0.05 


101.4±6.0 


21.0 


21.3 
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Table 3 — Continued 



No. 


name 


My 




log L/L Q 


R 


M spec 








mag 


mag 


dex 


Re 


M e 


M e 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 
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Table 4. Mass-metallicity Relationship of Galaxies 



Galaxy log M stars /MQ [Z] source 
(1) (2) (3) (4) 



M81 


10.93 


0.08 


a,b 


M31 


10.98 


0.04 




MW 


10.81 


0.00 


%,h 


M33 


9.55 


-0.15 


hj 


NGC300 


9.00 


-0.36 


k,l 


LMC 


9.19 


-0.36 


i./ii 


SMC 


8.67 


-0.65 


i,n,o 


NGC6822 


8.23 


-0.50 


hP 


NGC3109 


8.13 


-0.93 


i>« 


WLM 


7.67 


-0.87 


i,r 


Sex A 


7.43 


-1.00 


i,s 



de Blok et al l 120081) 
b this work 



IChemin et al . 2009) 



Przybilla et al. 2008b) 



e Trandle et al. (2002) 



Smartt et al. (2001) 



fSofue etal. (2009) 



' lPrzvbilla et al l (2008a) 



Woo et a l. (2008) 



Uetal 



(2009) 



Kent 



1987) 



Kudritzki et al. (2008) 
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n Hunter etal. 2007) 



Schiller 



2010) 



Trundle & Lennon 2005) 



1 Venn et al 



2001) 



c Evans et al 



2007) 



Urban eia et al. 



2008) 



Kauferetal. 2004) 



